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Abstract 

EuBq is a magnetic semiconductor in which defects introduce charge carriers into the conduction 
band with the Fermi energy varying with temperature and magnetic field. We present experimental 
and theoretical work on the electronic magnetotransport in single-crystalline EuBg. Magnetization, 
magnetoresistance and Hall effect data were recorded at temperatures between 2 and 300 K and 
in magnetic fields up to 5.5 T. The negative magnetoresistance is well reproduced by a model in 
which the spin disorder scattering is reduced by the applied magnetic field. The Hall effect can 
be separated into an ordinary and an anomalous part. At 20 K the latter accounts for half of 
the observed Hall voltage, and its importance decreases rapidly with increasing temperature. As 
for Gd and its compounds, where the rare-earth ion adopts the same Hund's rule ground state as 
Eu'^^ in EuBg, the standard antisymmetric scattering mechanisms underestimate the size of this 
contribution by several orders of magnitude, while reproducing its shape almost perfectly. Well 
below the bulk ferromagnetic ordering at Tc = 12.5 K, a two-band model successfully describes 
the magnetotransport. Our description is consistent with published de Haas van Alphen, optical 
reflectivity, angular-resolved photoemission, and soft X-ray emission as well as absorption data, 
but requires a new interpretation for the gap feature deduced from the latter two experiments. 
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I. INTRODUCTION 



The binary compound EuBg crystallizes in a simple cubic lattice, with divalent Eu ions in 
their ^87/2 Hund's rule ground state at the corners of the unit cell and Bg-octahedra centered 
at the body-centered positions. With decreasing temperature, it orders ferromagneticallyi 
via two consecutive phase transitions at ~ 15.5 K and ~ 12.5 K, respectively^, the first 
of which has recently been interpreted as a phase separation between small magnetically 
ordered regions with mobile charge carriers and large disordered regions with localized mag- 
netic polarons^. A spin-polarized electronic-structure calculation in the local spin-density 
approximation for exchange and correlation (LSDA) correctly reproduces the lattice con- 
stant, the internal coordinates of the boron atoms and the size of the magnetic moment in the 
stoichiometric compound^. It also predicts the system to be a semimetal, with overlapping 
conduction and valence bands around the X-point of the Brillouin zone (BZ), in contradic- 
tion with the results of a combined study based on angle-resolved photoemission (ARPES) 
and bulk-sensitive soft X-ray emission (SXE) and absorption (XAS) spectroscopies^, which 
suggest a gap of at least 1 eV between the two bands. While, given the approximation used 
in the calculation, this discrepancy is not surprising, the fact that conduction band states 
were observed at all in the photoemission experiment, illustrates the fundamental problem 
one is faced with when trying to describe the transport properties of this system and of 
the hexaborides in general, namely that their behaviour is to a large extent determined by 
defects and impurities^. 

In this paper, we offer a quantitative analysis of magnetization, magnetoresistance and 
Hall effect data, obtained on one and the same sample, using a plausible model for the 
origin of the mobile charge carriers and a consistent description of the dependence of their 
concentration and their scattering rate on the applied magnetic field and on temperature. 
Explicitely, we consider EuBg to be a strongly compensated n-type magnetic semiconductor. 
Due to the merging of defect states, i.e., boron vacancy levels, with the conduction band, 
the latter acquires a certain concentration of charge carriers, as evidenced by ARPES^. In 
the paramagnetic state and in the absence of external magnetic fields, these are equally 
distributed over six pockets (three for each spin direction) centered at the X-points of the 
BZ. Due to the thermal ionization of deep trap states in the gap, the occupation of these 
states increases slightly with temperature. The exchange coupling between the conduction 
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electrons and the localized (spin-) magnetic moments of the Eu-ions leads to a lowering of the 
conduction band edge above Tc and to a splitting of the spin-up and spin-down bands below 
T(r^. It is also responsible for the so-called spin-disorder resistivity^iiSiiLiiii^ which, in a 
semiconductor, is strongly dependent on the degree of spin-polarisation of the conduction 
electrons and on the concommitant redistribution of the charge carriers between the spin- 
up and spin-down bands-. We have attempted to model our data on the anomalous Hall 
effect with the mechanism suggested by Kondo^^ for the case of gadolinium metal, where the 
trivalent Gd ions adopt the same 4f'' configuration as Eu^"*" in EuBg. We have generalized 
this approach to include inelastic spin-flip processes. As in this previous work, we find 
that the shape of the calculated field- and temperature dependence of the anomalous Hall 
resistivity curves matches that of the measured ones almost perfectly, but the magnitude of 
the effect resulting from the calculation is several orders of magnitude too small, if reasonable 
values for the parameters are used in the theory. An alternative mechanism proposed by 
Maranzana^^, namely the interaction between the orbital motion of the conduction electrons 
and the localized moments, leads to the same functional dependence of pn on temperature 
and applied field, but with an even smaller amplitude. 

II. SAMPLE AND EXPERIMENTAL SETUP 

The single-crystalline sample of EuBg was prepared by solution growth from Al flux. 
All measurements were made using the same platelet type specimen with dimensions of 
approximately 4.8 x 5.2 x 0.25 mm^. The room temperature lattice constant of 4.185 A 
was evaluated from X-ray powder diffraction data, using a least-square refinement based on 
Cohen's method, with the software Xlat^^. A Si spectrum served as the internal standard. 
Gold wires with 25 /im diameter were contacted to the sample with silver epoxy. All volt- 
ages were measured with a four-probe, low-frequency ac technique in the ohmic regime. The 
transverse magnetoresistance, which in the following is always referred to as the magnetore- 
sistance, and the Hall voltage Vh were measured in a configuration where the external field 
Bat between and 5.5 T, was oriented perpendicularly to both the applied current and the 
measured voltages, thus orthogonal to the platelet. The extended temperature range was 
covered by using a conventional '^He cryostat. Magnetization measurements were made in the 
same geometry with a commercial superconducting quantum interference device (SQUID) 



3 



magnetometer, reaching temperatures between 2 and 330 K and magnetic fields up to 5.5 
T. 



III. ELECTRICAL TRANSPORT ABOVE 20 K 
A. a.) Theory 

The unperturbed conduction electron levels are approximated by parabolic bands, cen- 
tered at the X-points of the BZ, with spin-independent and, for convenience, isotropic effec- 
tive masses m*. The energies of the band-bottoms are specified by eg, where s = ± labels 
the spin of the electron. The corresponding Bloch states are denoted as |i k s) , where i = 
1, 2, 3 specifies the X-point from which k is measured. 

Following Haas^, we denote the eigenstates of the system of magnetic Eu-ions by \a) 
and their occupation probabilities by Wa, so that, for example, the equilibrium value of the 
z-component of the spin of the ion located at the site i?„ is given by 



Our magnetization measurements show that above ~ 20 K and up to fields of 6 T, this 
quantity is well described by molecular field theory and in the following, we shall assume 
that in this temperature range, this also applies to higher-order correlation functions. 

The exchange interaction between the magnetic moments and the mobile charge carriers 
has the form 



where the sum is over all unit cells in the crystal, and s is the spin of the conduction 
electron at V. The range of is determined by the radius of the 4f-shell. 

To first order in perturbation theory, this interaction produces the following modification 
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in the energy eigenvalues of the band electrons 



Ae^^-'( k ) = '^^Wa{i k s; a\Hi\i k s; a) 



^ ^ _ ^ - 

= — ^^(^ ^ \JC^ ~ Rn)\i k) ■ (s|'s*|s) ■ Wa{a\Sn\ce) 

n=l a 

For a collinear ferromagnet, such as EuBg, the average value embodied in the last sum 
points along the magnetization axis, which we choose as the quantization axis for the con- 
duction electron spins. Furthermore, all magnetic ions being equivalent, we can define 

X^ = N{z^\J{V)0)^J , (4) 

where, in the last step, we have used the effective mass approximation. 
The spin-dependent energy shift then takes the form 

Ae« = -i.Ja , (5) 

where the mean ionic spin a=S{M/Msat) depends on temperature and magnetic field. 
Msat=— (N/V)g/XBS is the saturation magnetization per unit volume, V is the volume of the 
sample, fiB is the Bohr magneton, and g is the g-factor of the magnetic ion with spin S. 

The second order correction is independent of the spin of the conduction electrons in the 
absence of a net magnetization and is completely dominated by the first order splitting (see 
eq. (0)) otherwise, so that it may safely be neglected in the discussion of magnetotransport 
properties. 

The generic form for the matrix element of the interaction 7-^2 responsible for the anti- 
symmetric scattering between states with wave vectors k and k' isf^^ 

{V±-a\n2\k±-a) = iC (k X fc') ■ z ^ e*(^-^)^(a|5„,|a) , (6) 

n=l 

where C ( ^-r^Tt 1 is an even function of the scattering angle 6'-+->. This describes the scat- 

\ k k J k k 

tering along the direction perpendicular to k and to the magnetization for spin-conserving 
transitions, and does not affect the position of the energy bands to first order in perturbation 
theory. 
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For the coupling between the orbital motion of the conduction electrons and the localized 
spins we write, as proposed by Maranzana,— 



n=l 



and 



\ k k J m*V C0S^6-^-rt 

k k 

For the mechanism suggested by Kondoli to explain the anomalous Hall effect in gadolin- 
ium, which involves a virtual excited state with one electron less in the 4f shell, 

^(M = 2:^bv ' 

where A is the spin-orbit radial integral for 4-f electrons in the 4f^ configuration, Vi is the 
mixing matrix element between the / = 1 component of the plane wave factor e^^ of the 
Bloch function (which is modulated by a function m-^ with the periodicity of the lattice and 
of pure d-character around each Eu-site) and a 4f-orbital, and A_ is the minimum energy 
necessary to excite one electron from the 4f'' configuration to the Fermi level. 

The non-periodic part (Ti'i) of Tii, obtained by replacing the z-component of Sn, Snz, by 
Snz — cr in eq. Q, and H2 induce transitions between different Bloch states, the probability 
of which is given to lowest order by Fermi's golden rule as 

P^^\ts-a\Vs'-a') = '^6{e?-.-e?:)\{Vs'-,a'\n[+n2\'ks-,a)\^ . (10) 

The energy transfer between the conduction electrons and the spin-system has been 
neglected (quasielastic or quasistatic approximation), which is justified as long as the typical 
excitation energy of the latter is smaller than the thermal energy^, as is certainly the case 
above Tc. The rhs of eq. (jl(J|) contains 4 terms. The two cross products cancel due to the 
fact that the matrix elements of 7^2 are imaginary. From the definitions (jH)) and it is 
easy to see that the contribution from the square of the matrix element of 7^2 is a priori 
irrelevant in the case of Maranzana's mechanism and contributes less than one percent to 
the total transition probability ()10|) if a physically reasonable value is used for the ratio 
Vi/A^ in Kondo's model. Therefore we expect that the temperature and magnetic field 
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dependence of the resistivity is controlled entirely by the exchange interaction between the 
conduction electrons and the localized moments. 

Two types of transition have to be considered, those without spin-flip by 



(j:?±;a|7^;|zt±;a) = T^^Jg_,.eW^-^-?^+^^^^^^^ , (11a) 

n=l 

and those with spin-flip by 



{jVt; a ± l\n[0±; a) = -^f2 J'h ' ^'^^'^^^^^^(fi^ ± ^S^W) , (lib) 

n=l 

where we have introduced the spin-raising and lowering operators = Snx ± iSny at 
the site Rn, and 

jli^ = N-{jV\Jir)\t~k) . (lie) 

For intravalley transitions {i = j), the momentum transfer is small, and we can set Jl^^, 
equal to J defined in Eq. (0} above. The short range of the exchange integral in real space 
implies that the matrix element for intervalley scattering will not be much reduced with 
respect to J. Fortunately, the short wavelength of the associated spin fluctuations and, in 
particular, the small range of scattering angles available for this process, allow us to neglect 
it. The same argument can be used to dismiss intervalley scattering in 

To lowest order, the transition probabilities associated with the matrix elements ()lla|l 
and ()llb|) are then given by 

P(^)(t±; a\V±; a') = ^5(ei - e|) (^V ■ eW^-^)(^-^)] ((5.. - a)(5„,. - a)) , 

n,n' 

(12a) 

and 

p(^)(t±;a|^T;a') = ^5(4 - e|) (^X • eW^-^)(^-^)l . (12b) 
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Following Haas^, we express the spin correlation functions appearing in Eqs. ()12all2b|) 
in terms of the generalized susceptibility per unit volume 

xHt) = ^e[^7 (^-^)l X - {S^.){S^,)} , (13) 

n,m 

where i,j = x,y,z. For a simple cubic (lattice constant a) coUinear ferromagnet, with 
one magnetic atom per unit cell, x*"* is diagonal and can be written as 

X\t) = U,r' + Aq']-' (14) 

for small values of ~q, where Xh is the susceptibility per unit volume of a single domain 
in a homogeneous magnetic field and 

A = (151 

The transport relaxation rate for a Bloch state \ k ±), with energy is then given by^ 



Inserting the explicit form of the susceptibilities and performing the sum over k' , we 
finally obtain, for spin disorder scattering, 



1 



m*kBT ( J 



^ y k 



X 



NgfiB 
m*A 



) + ^4^^^ 



2 , / ' Ahv ^1 ^1 \l k k 



(17) 



'^''^ fc k V k k- 

In the molecular-field approximation, we have^ 



M = M,atzBs{gfiBS\F\/kBT) , 



(18) 
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where Bs is the Brillouin function for a spin 5* acted upon by the effective field 

F=fiojM + X . (19) 
Here 7 is the molecular field constant and Ba is the applied external field, so that 

xl = xl = M/Ba , (20) 

and 

Xl = [(M,,, (dBs/dF))-' - f^olV' . (21) 

The antisymmetric scattering responsible for the observed anomalous Hall effect has its 
origin in the matrix element given in eq. (jHI). For it to appear linearly in the transition 
probability, we need to compute the latter to third order in the matrix element o^^i^^i^^ , which 
leads to 



p(=^V"fc ±; a\V±; a) ^ ^5(ei - et] 



{k' ±; a\H2\ k ±; a) ■ 



( k ±; a\ni\k"s"; a"){k"s"] a"\ni\k' ±] a) 



et - eC + i6 

k k" 



(22) 



where the real part is derived by use of the identity 



1 



k" 



V 



1 



k k" 



(23) 



In contrast to Kondo and Maranzana, we allow for spin-flip exchange scattering to and 
from the summed-over intermediate states, but the inelasticity in energy has again been 
neglected. Inserting the appropriate matrix elements, we obtain 



P^^\'k±;a\k'±; 



27r / J 



a 



h \2N 



X 



± s" 
k k" 



Df{k,k\k") 



(24) 



k'\s" 
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with 



".p.? (25) 
■ [{SUSp, - a){S,, - a)) + {Sn.S;S^)] . 

In the spirit of molecular field theory, we now assume that the sum over three-spin 
correlation functions can be limited to those terms in which n = p or n = q, and then make 
the following decoupling 



^i(k"-k')-(Rp-Rg) _^ ^i{k"- k)iRp-Rq) 



p,q 



(26) 



• [{{S,^-a)iS,,~a)) + {S;S^)], 
which, with the definition of the generalized susceptibility (fT!?jl . leads to 



D'i\k,k',k") = Vj^^a 



X^(k"~k')+x'ik"- k)+2x%k"-k') + 2x%k"- k) 

(27) 

The sum over intermediate states in Eq. is now identical to the one appearing 

in the expression for the transport relaxation rate of Eq. p6|) . with = or for 
spin conserving or spin-fiip transitions, respectively. This allows us to write the transition 
probability for skew scattering in the compact form 

P(3)(t±; a\'P±; a) = 2C (o^p) 6{e^ - e|) ■ {k x k') ■ z ■ ■ a , (28) 



Tie: 

^ k ■ 



where we have used the fact that the relaxation rate of a particular Bloch state only 
depends on its energy. Given an interaction with the above angular dependence, we can use 
the exact result derived by Fert-'^^ for the Hall resistivity ph (his Eq. (15)), to define the 
relaxation rate for antisymmetric scattering 



V 



f d^k'P^^\~k±- a\y±- a)k' ■y = C- k—^ ■ a , (29) 
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where, for the interaction proposed by MaranzanaM 



= ^^f^ = 1.794- 10-^^ m . (30) 
10 



For the mechanism, suggested by Kondo^^ to explain the anomalous Hall effect in gadolin- 
ium metal, 

^ ~ Qir^S^lh? ' ^'^^^ 
where Vlceii is the volume of the unit cell. 

In our comparison with experiment, the constant C in eq. ()29|) will be treated as a free 
parameter. 

B. b.) Experimental Results and Analysis 

1. Temperature dependent electrical resistivity in zero magnetic field 

The temperature dependent electrical resistivity arises from several scattering processes, 
which we shall, as usual, consider as independent (Matthiessen's rule). From the theory 
developed above we can calculate the conductivity of the coupled spin-up and spin-down 
charge carriers in the presence of spin-disorder scattering only according to 

-UB^^T) = [jf^^^)^^ - 4)9Hef-^de (32) 

with the Fermi function 

/o(e) = ^ , (33) 

1 + exp{j^) 

where C, stands for the chemical potential. The density of states for each spin orientation 

is 

The relaxation times T^(e) = '^"(ei) ^^re given by Eq. ()17|1 . In the absence of a magnetic 
field and in the temperature range considered here, the system is unpolarized (ed" = eg), 
and the relaxation time is the same for both spin orientations. We can then define the 
spin-disorder resistivity as 



p,,(0, T) = (a+ (0, T) + a;,(0, T))-' . (35) 
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In order to calculate the contribution pph of the electron-phonon interaction to the re- 
sistivity, we use the model that was recently suggested by Mandrus and collaborators^^ for 
LaBg. The electrons are assumed to be scattered by localized low-frequency Einstein oscil- 
lators, corresponding to the almost independent motion of the rare-earth ions in their boron 
"cages", as well as by Debye-type phonons due to the collective motion of the boron frame- 
work. In a first step, we apply the model, described in detail in ref.— , to fit the resistivity 
data of YbBg^^. In this compound, the Yb cations also adopt a divalent configuration but 
they carry no magnetic moment. We then renormalize the obtained Einstein frequency by 
the square root of the mass ratio between Yb and Eu, leading to ^£;=168 K for EuBg. The 
Debye frequency (6'£)=1160 K), to which the results are not sensitive to start with, is left 
unchanged. Next, we have to account for the resistivity pd arising from the scattering of the 
conduction electrons at point defects. We anticipate the charge carrier density to be high 
enough to efficiently screen the latter and therefore, the corresponding relaxation rate can 
be considered as temperature-independent. The total resistivity is then given by 

p(0,T) =p,rf(0,r)+pp/,(0,T)+prf(0,r)+p,o„t , (36) 

where pcont is a (small) contribution arising from non-ideal electrical contacts to the 
sample, and which we assume to be independent of temperature and magnetic field. 

In the next step we compare eq. (jHSl) with the measured temperature-dependent re- 
sistivity. To begin with, we postulate that the mobile charge carriers in the conduction 
band originate from the transfer of electrons from doubly and singly occupied levels of 
Bg-vacancies. The existence of such defects has been invoked by Noack and Verhoeven^S to 
explain their gravimetric data on zone refined LaBg. Their formation energy has been shown 
to be substantially smaller than that of six widely separated B-vacancies^l. An excellent fit 
is obtained in the range 40 K < T < 100 K with a constant carrier concentration of 1.4-10^^ 
m~^ or 10~^ / unit cell, which corresponds to a Fermi energy Ep of 54 meV. At elevated 
temperatures, the experimental data suggest that electrons from a narrow "band" of defect 
states which, for reasons that are elucidated below, we associate with compensating ionized 
acceptors in the form of Eu-vacancies, start to populate the conduction band. The experi- 
mental data is well reproduced if we assume a concentration of 6-10^^ m~^ defect levels, with 
a lorentzian energy distribution centered at 19 meV below the conduction band edge (i.e., 
73 meV below Ep) and a full width at half maximum of 9 meV. Finally, our fit requires the 
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density of mobile charge carriers to increase by 40 % as the temperature is reduced from 
40 K to 22.5 K. This increase can be explained by an early onset of magnetic short range 
order—, nucleated by the presence of defects, which locally reduces the activation energy of 
the donor states. The different components of the resistivity, their sum, and the measured 
curve are displayed in Fig. ^ In contrast to earlier work^^, where the contribution from 
electron-phonon scattering to the room temperature (RT) resistivity was estimated to be 
less than 3 %, our analysis shows that this mechanism actually dominates above 125 K and 
is responsible for over 60 % of the total resistivity at room temperature. Due to the small 
size of pd + Pcont, an unambiguous estimate of the contact term is not possible at this stage 
and requires the analysis of the magnetoresistance given below. The variation of the charge 
carrier density with temperature is summarized in Fig. |21 

2. Magnetoresistance 

The magnetization of our sample as a function of the applied magnetic field is displayed 
for a large number of temperatures in Fig. El which also shows the results of a fit using 
Eq. (fTH|) to all measured temperatures and fields above 30 K. The latter yields a saturation 
magnetization of (8.83 ± 0.04) • 10^ A/m, in excellent agreement with the value of 8.86 ■ 10^ 
A/m expected for divalent europium, and an effective molecular field parameter 7 = 5.15 ± 
0.05. Besides the Weiss field, 7 contains the Lorentz field (7^ = (1/3)), negligible in higher 
Tc materials, and the demagnetizing field {pfo ~ —0.93 for our geometry). The Curie 
temperature of an infinite size bulk sample is determined by the first two terms and amounts 
to 13.6 K, close to the temperature at which neutron scattering experiments^i reveal the 
onset of spontaneous magnetic order. 

The parameters also allow to calculate the longitudinal and transverse susceptibilities 
using eqs. (j7T|) and (^01) . respectively, as well as the shift of the bottoms of the spin- up and 
spin-down conduction bands, induced by the non-zero magnetization, via eq. (0). The latter 
leads to a redistribution of charge carriers between the two bands which, in turn, requires 
an adjustment of the chemical potential with respect to the band minima. The Eu-vacancy 
levels will also, to a lesser extent, be affected by the magnetization. The spin-down states 
will rise in energy and progressively empty themselves into the (spin-up) conduction band. 
At some temperature-dependent value of the field, the latter will merge with the spin-up 
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Eu-vacancy states. Due to the Pauli principle, the transport properties will not be affected, 
however. 

The resistivity may now be calculated as follows. For all values of the applied field and 
temperature, we define two average relaxation rates 



1 1 n^e^ 



(37) 



due to spin-disorder scattering. According to the model of Mandrus et al.^^ the electron 
phonon relaxation rate is proportional to the Fermi velocity and hence, because vp ~ n^^'^ 
we can write 

. / ± X 1/3 T 2 / ± \ 1/3 



'^ph KVI^I/ 'ph 



Finally, we assume that the (weak) scattering by point defects is field- independent, which 
leads to the total average relaxation rate 



m* 



— + n[^] pp,(0, T) + np,(0, T) 



(39) 



and to the total resistivity in the presence of a magnetic field 



p{Ba, T) = ^ {n+T+ + n-T-) + p,ont (40) 

In passing we note a spin-polarization of the itinerant electrons, arising from a redistri- 
bution of the charge carriers between the spin-up and spin-down bands, resulting from the 
opposite shifts of the band edges and described by eq. (0). These shifts alone leave 
the density of mobile charge carriers constant but the negative shift of the majority band 
leads to a transfer of carriers from localized defect to itinerant band states. 

The free parameters in the model are the exchange coupling constant J, the effective 
mass m*, the contact resistivity pcont and the charge carrier density ntot{Ba). The best 
agreement with experiment is obtained for J = 0.18 eV, m* = 0.22 ■ rrie, where rrie is the 
free-electron mass, and Pcont = 1.5-10^'' Qm. A ten percent (correlated) variation of the 
parameters still produces reasonable results. Our optimum value for the exchange coupling 
constant is very close to the one quoted by Rys et al.® for divalent europium (0.188 eV) and 
our value for m* compares well with the density of states mass tudos = 0.26-me yielded by 
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the LSDA bandstructure calculation^. The absolute values of the carrier densities n±{Ba) 
depend strongly on J and m*, but their relative changes are identical for all parameters. 

Figure m displays the measured curves for p{Ba,T) at 22.5, 40, 60, 80, 125 and 175 K. 
The solid lines represent the calculations at the corresponding temperatures and fields. We 
note a perfect agreement at temperatures above 60 K. At 22.5 K and fields less than ~ 
2 T, strong polarization effects induce a substantial variation of n{Ba) which is difficult 
to model. Nevertheless, the calculated curve reproduces the measured results to within 5 
%. For stronger magnetic fields the measured curve for p{Ba) decreases more slowly with 
increasing magnetic field, reflecting a further reduction of the spin-disorder scattering and 
an increase of the charge carrier density. Eventually, p{Ba) flattens out and subsequently 
increases slightly towards the highest fields, which we interpret as the onset of conductivity 
through a second band, described in more detail in the section on the low-temperature 
transport. The charge carriers in this second band are holes with a concentration increasing 
from at 4 T to n,, = O.S-lO^^m'^ at 5.5 T. 

In Figures El a and b we display the electron densities ri^ and n~ in the spin- up and spin- 
down band, respectively, for the same temperatures and fields for which p{Ba) was calculated. 
With decreasing temperature, the polarization effects lead to a stronger enhancement of 
and a corresponding reduction of n~ with increasing field. At 60 K n' vanishes at ~ 4.5 T, 
leaving a fully polarized conduction band. At 40 K, n~ vanishes at 3 T and at 22.5 K already 
at 1.4 T. As mentioned before, an increasing concentration of holes has to be introduced 
below 40 K in order to explain the high field {Ba > 4 T) data. The only plausible mechanism 
for this to happen is that, e.g., at 22.5 K and 4 T, the top of the spin moment up valence 
band which, according to the calculation of Massidda et al.— should experience an (upward) 
shift of the order of 15 percent of that of the bottom of the conduction band, touches the 
Fermi level Ep. Note that this is not in contradiction with the existence of ionized (i.e., 
occupied by electrons) acceptor states below Ep. 

3. Hall effect 

In a magnetic conductor, the Hall resistance consists of two contributions, namely the 
ordinary part p^'^, due to the Lorentz force ev' x B acting on the electrons, and the 
anomalous part p^"^, which results from the antisymmetric scattering of itinerant charge 
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carriers by the disordered local moments on the Eu-ions^^. The spin-flip exchange scattering 
mixes the states of the spin-up and spin-down conduction bands, which can therefore be 
considered as a single entity. Therefore the ordinary Hall resistivity is related to the total 
density of mobile charge carriers, ritot by the usual relation 

p'S' = -^^-B , (41) 

with 'B =bI + /io(l - 7d)M. 

In Fig. Elwe display the measured Hall resistivity pn as a function of applied field for 22.5, 
60 and 125 K, together with p^'* computed with the charge carrier densities obtained from 
the fit to the magnetoresistivity data. The difference between the measured Hall resistivity 
and p'S'^ is largest at low temperatures and in small fields, where the spin up and spin down 
carrier densities are most sensitive to fluctuations in the magnetization (see Fig. E)). 

We attempted to model this difference, which we interpret as the anomalous Hall resis- 
tivity p^'*^, using the relaxation rate for antisymmetric scattering given by equation (|29|1. 
Treating (7 as a free parameter in a fit to p^"^ with p^°^ = {a^°'^'~^ + a^"'^'~)~^ , and inserting 
the Hall conductivities obtained from eq. (jH^ with Tsd replaced by Tas, yields the curves 
shown in Fig. [7|and the optimum value C ~ (6.6 ± 0.5)-10~^^ m, exceeding that of C*^ (eq. 
(j30|) ) by more than four orders of magnitude. Using the value of m* obtained from our fit to 
the magnetoresistivity, the spin-orbit coupling constant A = 164 meV for the intermediate 
4f^ configuration^^ and the lattice constant a = 4.185 A, we can write the corresponding 
coefficient for the Kondo mechanism as 

= 1.77-10-1=^^ m . (42) 

For to adopt the optimum value obtained from the fit to pn — P°h'^ would require 
the ratio Vi/A_ to be of the order of 20, which is utterly unrealistic. According to X-ray 
photoemission experimentsi^i, the lower bound on A_ (''Fq final state) is ~ 0.7 eV and Vi is 
expected to be smaller. It appears that discrepancies of that order are the rule for systems 
with half-filled 4f-shells in their ground state, such as gadohnium and Gd compound a^^i^^ . 



16 



IV. LOW-TEMPERATURE TRANSPORT 



At temperatures below 8 K, electron-phonon scattering is negligibly small, and the only 
contributions to the resistivity are p^d, Pd and pcont- In figures (jHI) and 0, the data for 
the magnetoresistance and the Hall resistance are plotted for 2, 4 and 8 K. At 2 K, p{Ba) 
increases by a factor of ~ 7 between and 5.5 T. In fields exceeding 1.5 T, p{Ba) is nearly 
quadratic in Ba for all three temperatures. This observation strongly suggests that two 
bands with oppositely charged carriers participate in the conduction of electrical current. 
For Ba > 1.5 T, we therefore use a standard two-band model^ to simultaneously describe 
Ph^Bo) and p{Ba). This leads to 

_ Ripl + -R2P1 + RiR2{Ri + R2)Bl 

{p, + p,Y + {R, + R.fBl ^^'^^ 

and 

_ PlP2{Pl + P2) + {piRl + P2Rl)Bl , 

p - {p, + p2f + {R, + R2yBi + ' 

where i?2, p2 and p2 depend on temperature and on the applied field. The Hall 
"constant" R is the proportionality factor between pniBa) and Ba'-, Ri and R2 are the Hall 
"constants" for the conduction and the valence band, respectively. 

At 8 K and in zero external field, the ordered Eu moment is of equal magnitude^^ as 
the net moment per Eu ion in the field direction at 22.5 K and Ba = 5.5 T. Hence we 
expect to find the same concentration of electrons rif, ^ 6.1-10^^ m~'^ in the (fully polarized) 
conduction band and the same value of Ri{8 K, T) ^ — 1.0-10"^ m^A~^s~^ in both cases. 
Similarly Uh ^ 0.8-102^ m'^ and i?2(8 K, T) ^ 7.7-10"'^ m^A^^s'^ We determine pi and 
P2 under the same conditions as follows. First we note (Fig. (jSJ) that the contribution from 
spin-disorder scattering to the total resistivity is negligibly small for 5^ > 1.5 T. For applied 
fields in excess of this value pi (P2) is therefore entirely due to the scattering of electrons 
(holes) by point defects, and is proportional to the density of electrons (holes), with no 
explicit dependence on Ba- This allows us to extrapolate this contribution, which we call 
Pid (P2d), to zero applied field as follows 



n (22 5K OT) 

Pu{8K, OT) = pid(22.5i^, 5.5T) = ^ ^^2 5K 5 5T) ' '"i'^^22.5ir, OT) ^ 0.66 ■ lO'^fim (44) 
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where the carrier densities can be read off Fig. (0), and pid{22.5 K, T) is obtained from 
Fig. (P) and pcont determined in the previous section. From the two band model and in the 
absence of spin-disorder scattering it follows that P2d{8 K, T) 1.44- 10~^ Qm. Considering 
the ratio of the effective masses and the carrier concentrations for the two bands, we find 
that the relaxation time of the holes is approximately three times shorter than that of the 
electrons. Enhancing the applied field from to 5.5 T induces a monotonous enhancement 
of the ordered Eu moment and thus the magnetization. This in turn enhances the overlap 
between the valence band and the conduction as well as the donor spin-up bands, leading 
to a net increase in the density of mobile charge carriers. From the fit of our experimental 
data to eqs. (HHafl and (l43Hll for Ba > 1.5 T we obtain = G.T-IO^^ m'^ and Uh = 6.1-10^^ 
m~'^ at magnetic saturation. The growth rate is roughly proportional to {M / Mgat)^^'^ , as 
expected for parabolic bands. The residual spin disorder resistivity in zero field amounts to 
less than 1-10"^ fim (see Fig. (jH))). 

Below 5 K the elementary excitations of the system of magnetic Eu ions are spin-waves^i. 
From a comparison of the zero-field resistivities at 2 and 4 K we see that the scattering of 
the charge carriers by these collective modes, which should be proportional to T^, can be 
neglected. The field dependence of the (Hall) resistivity at these two temperatures is again 
well reproduced by the two-band model, on which we have imposed the constraint that the 
resulting values for Ri^2 and pi.2 at full magnetization are the same as at 8 K. 

V. DISCUSSION 

In this paper we offer a consistent, quantitative description of the magnetoresistance and 
the Hall effect in EuBg over a wide range of temperatures above and below the magnetic 
phase transition. From our analysis, the following picture of the electronic structure of this 
compound emerges. For T > 20 K and no applied magnetic field, EuBg is a heavily (self-) 
doped, strongly compensated n-type semiconductor. The donors are Bg vacancies, whose 
energy levels form a narrow "band" centered above the chemical potential, and possibly 
trivalent impurities, which we have neglected. The acceptors are cation vacancies, always 
present in the hexaborides, whose energy levels also form a narrow "band" just above the 
top of the valence band (Fig. ^Ja). At 22. 5K, the intrinsic band gap is of the order of 
14 meV and the chemical potential lies ~ 66 meV above the bottom of the conduction 
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band. Upon application of a magnetic field, the Eu-ions acquire a net moment in the field 
direction. The latter couples to the conduction electrons through the exchange Hamiltonian 
Til, leading to a splitting of the conduction band proportional to the magnetization in first 
order. The acceptor levels, the wave functions of which are mainly composed of (1^2 orbitals 
reaching out from the neighboring cations, suffer a splitting of similar size, while that of the 
donor levels and especially that of the valence band, which both couple to the Eu moments 
only indirectly through hybridization with the cation's d-orbitals, will be smaller—. These 
splittings lead to a redistribution of the charge carriers between the different bands and, 
in our model, at a temperature dependent critical value of the applied field, e.g., 1.4 T at 
22.5 K, the spin-down conduction band has emptied itself completely. At that point, the 
magnetization has reached only 25 percent of its saturation value. As the applied field and 
the magnetization are further enhanced, the top of the spin-down valence band moves closer 
to the selfconsistently determined chemical potential (Fig. (fTUlb)). until it crosses it. At 5.5 
T {M/Msat ~ 0.65) the situation is that of Fig. (fTUl c). Once the saturation magnetization 
has been reached, the carrier densities in the valence and the conduction band stay constant. 

This picture is consistent with the observation of two ellipsoidal pockets in de Haas van 
Alphen (dHvA) and Shubnikov de Haas experiments, performed at fields above 5 T-^*^. 
The electrons and holes being in different spin states, their Bloch functions cannot mix. 
It also offers a natural explanation for the weak temperature dependence of the dHvA 
frequencies, even accross T^^, since these are only affected by deviations of the 4f-electron 
based magnetization from its saturation value. In view of the sensitivity of the system to 
defects, our carrier concentrations are in very reasonable agreement with the ones quoted in 
refs.— ancl^. 

A further test of our interpretation is provided by the reflectivity experiments of Degiorgi 
and collaborators^*^. In Fig. ^2 we display the bare plasma frequency 



as a function of temperature in zero field, calculated with the carrier densities obtained from 
our fits and the optical masses provided by the band structure calculations of ref;^ (m^^* = 
0.24me, nf^^ = 0.29me). Our results compare well with the data of ref.— . In particular, we 
reproduce, even quantitatively, the steep rise of Up below Tc- The decrease of the computed 
plasma frequency at temperatures between 300 K and 100 K is consistent with the observed 




(45) 
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red shift of the plasma edge in ref.— . 

Fig. ^1 shows the dependence of the plasma frequency on magnetization obtained from 
our model. In ref.—, uj"^ was found to be proportional to M for 1.6 K < T < 35 K and 
T < -Ba < 7 T. Again, our results are compatible with this behavior, but suggest that the 
relation between the two quantities may be more complex. 

Finally we compare our model with the ARPES and bulk-sensitive XAS and SXE data 
of ref.—, which were obtained in the temperature range between 20 and 30 K. The existence 
of an X-point electron pocket was assumed a priori in our theoretical ansatz, relying on the 
validity of these experimental results. The authors of ref.^ attribute the feature labelled 
"band 1" in their paper to the emission from the valence band. The fact that its dispersion 
is much weaker in EuBg than in CaBg and SrBg, and the value of its binding energy of ~ 1.2 
eV, lead us to interpret it as an emission from the Eu 4f-shell. According to our model, the 
emission from the valence band should start at a binding energy of ~ 0.1 eV, which is not 
seen in ARPES. This may be due to the fact that the exposed [100] face consists of metal 
atoms only and that the electrons, originating from the boron network, cannot escape from 
the solid at the given photon energies. Another complication is the observed time-dependent 
surface relaxation^. The SXE and XAS data are consistent with our interpretation of the 
ARPES data. 

VI. CONCLUSION 

Although "real" EuBg is a heavily doped, strongly compensated, and therefore very 
disordered magnetic semiconductor, many of its properties can be satisfactorily described by 
a relatively simple model, taking into account the two main intrinsic sources of imperfections, 
namely Eu and Bg vacancies. Our microscopic treatment above Tc is limited to the range 
of temperatures where mean field theory can be used to describe the magnetic properties of 
the compound. We expect the behavior around Tc (and the critical temperature itself) to 
be sample dependent, as the (RKKY) coupling between magnetic ions is mediated by the 
conduction electrons whose concentration is a function of the magnetization. 

What is still missing is a plausible mechanism leading to the correct order of magnitude 
for the anomalous Hall resistivity in this and the other compounds formed by rare-earth ions 
with a half-filled 4f shell. In our opinion the order of magnitude discrepancy between theory 
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and experiment with respect to the anomalous Hall effect is not caused by underestimating 
the mixing matrix element Vi, as this would reflect itself in the resistivity as well. It must 
therefore be connected with the multiplicity of the intermediate states which can be coupled 
by the spin-orbit interaction. The lowest energy term for the 4f^ configuration is character- 
ized by L = 3 and S = 3, with a degeneracy of 49 in the absence of spin-orbit coupling. 
Whereas the unit operator, relevant for the resistivity, only has diagonal matrix elements, 
the spin-orbit operator, which splits the term into seven multiplet levels J = 0, 1,... 6, 
has matrix elements between states satisfying the selection rule AJ^ = 0, ±1, within every 
subspace corresponding to a given value of J. This leads to 133 possible transitions instead 
of 49, still not enough to account for the observed difference. We conjecture that the small 
hybridization between the europium (gadolinium) 4f orbitals and the boron sp^*2^ (gadolin- 
ium 5p^^) orbitals, suggested but overestimated by band structure calculations, opens the 
necessary extra channels. We hope that our results, which confirm previously established 
discrepancies, will encourage more theoretical work on this long standing problem. 
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FIG. 1: The total measured resistivity p of EuBg is represented by open circles. The calculated 
contributions due to scattering by phonons and magnetic excitations are shown by open squares 
and closed triangles, respectively. The closed diamonds represent the combined contribution to the 
resistivity of the scattering by point defects and of the non-ideal contacts. The solid line represents 
the total calculated resistivity in zero external field above 20 K. 
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FIG. 2: Itinerant carrier density n(T) in EuBg at high temperatures, obtained from the separation 
of the total resistivity into a magnetic, a phononic and an impurity contribution, (see fig. 1) 
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FIG. 3: Bulk magnetization M of EuBg as a function of applied magnetic field Ba oriented per- 
pendicularly to the platelet-shaped sample. All data for temperatures above 30 K, where fitted 
according to eq. (jlHj) . yielding the parameters Mgat and 7. The solid lines represent the mean 
field calculations for all temperatures using these parameters. Good agreement between this type 
of calculation and experiment prevails to even lower temperatures. 
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FIG. 4: Magnetorcsistivity of EuBe at 22.5, 40, 60, 80, 125 and 175 K, between and 5.5 T. The 
solid lines are the results of the resistivity calculations described in the text. 
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FIG. 5: The charge carrier concentrations n^{Ba) (a) and n^{Ba) (b) in the spin moment up and 
spin moment down band, respectively, at 22.5, 40, 60, 80, 125 and 175 K between and 5.5 T. 
Note the different scales of the y-axis. 
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FIG. 6: Measured Hall resistivities PHall{Ba) and calculated ordinary Hall resistivities pjj of EuBe 
at 22.5, 60 and 125 K between and 5.5 T. The empty symbols show PHall{Ba), the full symbols 
display p^"^. 
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FIG. 8: The measured (m) magnetoresistivity data p{Ba) for EuBg are displayed for 2, 4 and 8 K 
between and 5.5 T. The calculated (c) curves are obtained for the corresponding temperatures, 
using the two-band model captured by eqs. ()43a() and (|43b|) . 
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FIG. 9: The measured (m) Hall resistivity data puiBa) for EuBg are displayed for 2, 4 and 8 
K between and 5.5 T. The calculated (c) curves are again obtained, using the two-band model 
captured by eqs. (|43ajl and (|43b|) . 
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FIG. 10: (Color in online edition) Schematic electronic excitation spectrum of EuBg around the 
chemical potential fj, at T = 22.5 K. A possible arrangement of the conduction band, the valence 
band, an acceptor and a donor defect band is plotted for a: Ba = T, b: = 4 T, c: Ba = 5.5 
T. From left to right in each panel: valence band, acceptor levels (cation vacancies), conduction 
band, donor levels (Bg vacancies). All energies are given in meV. The up-arrows and down-arrows 
denote the spin moment up and spin moment down subbands, respectively, ^(0) is the chemical 
potential for Ba = T, whereas /i denotes the chemical potential at the corresponding fields. The 
distribution of the charge carriers over the 4 different bands is explained in the text. Note that in 
b) and c) the conduction subband for the down moment lies far above the chemical potential and 
is thus irrelevant for our purposes. 
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FIG. 11: Bare plasma frequency tOp of EuBg, obtained from the calculated itinerant charge carrier 
densities, vs. temperature in zero magnetic field. Up is calculated using eq. (|45|) and the optical 
massed provided by the band structure calculations {m'/^ = 0.24me, 'rn'^^ = 0.29me). 
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FIG. 12: Squared bare plasma frequency ujp of EuBg versus the relative bulk magnetization M/Msat 
at 2, 8, 22.5, 40, 60 and 80 K. ujp is calculated using eq. and the optical massed provided by 



the band structure calculations {m"^^ = 0.24me, rnyj'^ = 0.29me). 
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